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ABSTRACT 


A  2-D  digital  filter  with  a  rational  frequency  response 
can  be  expanded  into  an  infinite  sequence  of  filtering  opera¬ 
tions.  Each  filtering  operation  can  be  implemented  by  con¬ 
volution  with  a  low-order  2-D  finite-extent  impulse  response. 
This  sequence  of  filtering  operations  can  be  viewed  as  an 
iteration  where  a  new  estimate  of  the  output  signal  is  formed 
from  the  previous  estimate.  If  a  convergence  constraint  is 
satisfied,  the  sequence  of  estimates  will  approach  the  desired 
output  signal.  In  practice,  the  number  of  iterations  is  finite. 
Consequently,  the  frequency  response  that  is  actually  realized 
by  the  iterative  implementation  is  an  approximation  to  the 
desired  rational  frequency  response. 
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A.  INTRODUCTION 


For  digital  signal  processing  tasks  that  are  not  per¬ 
formed  in  real  time,  a  processor  may  have  access,  at  least 
conceptually,  to  the  entire  input  signal  and  the  entire 
output  signal.  This  seems  to  be  particularly  true  for  image 
processing  tasks.  An  entire  image,  residing  somewhere  in  the 
computer  system,  is  used  as  the  input  to  a  processing  program 
that  generates  an  entire  output  image.  In  these  cases,  it 
seems  plausible  to  consider  examining  the  output  image  and 
using  it  in  conjunction  with  the  input  image  to  form  a  "better" 
output  image. 

The  iterative  implementation  for  2-D  digital  filters 
described  below  uses  this  familiar  concept  of  feedback  in  an 
interesting  way.  A  2-D  rational  frequency  response  can  be 
used  to  form  an  iterative  computation  involving  only  finite- 
extent  impulse  response  (FIR)  filtering  operations.  The 
rational  frequency  response  can  be  theoretically  implemented 
by  applying  the  iterative  computation  an  infinite  number  of 
times,  providing  a  convergence  criterion  is  met.  In  practice, 
of  course,  the  iterative  computation  is  applied  only  a  finite 
number  of  times,  so  the  frequency  response  that  is  actually 
realized  is  an  approximation  to  the  original  rational 
frequency  response. 

The  iterative  implementation  seems  well-suited  to  digital 
processors  that  have  the  ability  to  convolve  a  2-D  signal  with 
a  filter  kernel  of  limited  extent.  With  repeated  applications 
of  the  iterative  computation,  these  processors  can  be  used  to 
implement  kernels  of  much  larger  extent  and  greater  filtering 
ability. 
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Because  the  iterative  computation  involves  only  FIR 
filtering  operations,  the  iterative  implementation  can  be  used 
to  approximate  some  real  symmetric  2-D  rational  frequency  responses 
without  the  necessity  of  2-D  spectral  factorization  [1,2]. 
Consequently,  2-D  filter  design  methods  that  yield  real  rational 
frequency  responses  [1,2]  could  potentially  be  used  to  design 
filters  for  the  iterative  implementation.  Furthermore,  2-D 
real  rational  frequency  responses  could  be  obtained  by  applying 
a  McClellan- type  transformation  [3,4]  to  1-D  real  rational 
frequency  responses. 
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B.  DERIVATION  OF  THE  ITERATIVE  COMPUTATION 


A  2-D  rational  frequency  response  can  be  written 

H(w1,a>2)  =  A  ( ,  <»>2 ) /B  (1) 

where  A(w^,w2)  and  B(uj^,<i)2)  are  trigonometric  polynominals 
given  by 

A(u)1,(jj2)  =  E  E  a(m1,m2) 
ml  m2 

B(o)^/0)2)  =  E  E  b(n^,n2) 
nl  n2 

It  is  assumed  that  the  above  sums  have  a  finite  number  of  terms; 
the  arrays  afm^n^)  and  b(nlfn2)  have  finite  extent.  Without 
losing  any  generality,  we  can  also  assume  that  the  ratio 
A (to^ , u)2 ) /B , uj2 )  is  normalized  so  that  b(0,0)=  1. 

Now  we  define  the  function 

C(w^,u)2)  =  1-B(a)^,a)2) 

=  E  E  c(n1,n2)  exp[- jwj^-ju^n^  (4) 

nl  n2 


expl-j(D1m1- ju)2m2]  (2) 


exp  [- jd)1n^- j(jJ2n2]  (3) 
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where 

c(n1,n2)  =  -b(n1,n2)  for  (n1 ,n2 ) ^  (0 , 0)  (5) 

=  0  for  (n1 f n2)=  (0, 0) 

Using  the  definition  of  C  (aj^,a>2)  ,  we  can  write 

H(w1,o)2)  =  A(a)1,w2)/[1-C  (to1,o)2)  J  (6) 

If  we  define  X  ,  tn2 )  as  the  spectrum  of  the  filter's  input 
signal  x(n^,n2)  and  Y(m^,(o2)  as  the  spectrum  of  the  filter's 
output  signal  y(n^,n2)  and  we  assume  that  both  J  Y  f  ,  a>2 )  j  and 
|x(oj1,u>2)|  are  bounded,  then 

Y  (o)1,w2)  =  H  w2>X  (cj^  ,a)2) 

=  A(u1,a)2)X(w1,a)2)/tl-C(w1,uj2)  ]  (7) 

or 


=  A(to1,a>2)X(a)1,a)2)+C(a)1,co2)Y(a)1,to2)  (8) 

This  implicit  formula  for  the  output  spectrum  suggests  the 
iterative  formulation 

Y_i  =  0 

Yi(u>1,u)2)  =  A(a)1,w2)X(a)1,a)2)+C(uirw2)Yi_1(a>1,u2)  (9) 
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After  I  iterations,  we  will  have 


I 

Yj  (<i)i ,  u>2 )  =  A(oi^  ,U32)  X(u)^,io2)  E  C  (w^ ,  ) 

i=0 

Letting  I  goes  to  infinity,  we  see  that 

0° 

Yffi(u)1,u)2)  =  A(o)1,oj2)  X(o)1,a)2)  E  C1  (u^  ) 

i=0 

=  A(a)1,u)2)  X(a)1,a)2)/[l“C(u)1,a)  )  ]  (10) 

The  iterative  solution  will  converge  to  the  desired  filter 
output  if  the  geometric  series  in  equation  (10)  converges. 

To  ensure  this,  we  impose  the  restriction 


|  C  ((s>1  ^ u>2 )  I  <  1  (ID 

Because  and  c(nlfn2)  are  finite-extent 

arrays,  the  filtering  operations  necessary  to  implement  the 
iteration  (9)  are  FIR  operations.  Consequently,  the  rational 
frequency  response  (6)  can  be  implemented  exactly  by  an  in¬ 
finite  number  of  FIR  filtering  operations. 

If  condition  (11)  is  indeed  true,  then  it  is  possible 
to  show  that  Yi(co1,w2)  converges  uniformly  to  Yfu^,^)  as 
i  goes  to  infinity.  In  addition,  it  is  also  possible  to  show 
that  yi(n1,n2),  the  inverse  Fourier  transform  of  Y^ ,w2) , 
converges  uniformly  to  y(n1,n2)  as  i  goes  to  infinity. 


The  use  of  the  iterative  computation  (9)  can  be  visual¬ 
ized  as  a  simple  digital  filter  that  processes  a  sequence  of 
images  (or  vectors)  rather  than  a  sequence  of  numbers. 

Figure  1  shows  the  iterative  computation  as  a  first-order  feed¬ 
back  loop,  with  x^(n^,n2)  ^en°ting  the  input  sequence,  y^(n^,n2) 
denoting  the  output  sequence,  and  i  representing  the  iteration 
number.  The  "STORE"  operator  in  Figure  1  stores  the  result 
of  the  previous  iteration;  it  is  analogous  to  the  1-D  shift 
operator  usually  represented  by  z  To  obtain  the  output  signal 

y  (nl  ,n2  >  whose  Fourier  transform  YU^oj  )  is  given  by  equa¬ 
tion  (7),  the  input  sequence  is  taken  to  be  a  step 

function  in  the  iteration  variable  i: 

xi (ni ,n2 )  =  x(n1,n2)  for  i  >  0  (12) 

=  0  for  i  <  0 

As  suggested  by  equation  (9) ,  the  initial  condition 

Y-i (ni,n2)  =  0  (13) 


is  used  to  begin  the  iteration  at  i=0.  Then  the  desired  output 
y (ni ,n2 )  is  realized  as  the  steady  state  output  of  the  filter 
in  Figure  1. 

C.  STABILITY  AND  CONVERGENCE 


Stability  of  the  filter  H(o)1,w2)  and  the  satisfaction  of 
the  convergence  criterion  by  C(u)1,co2)  are  not  equivalent. 
Consider  the  filter 

H(w1,u)2)  =  l/(l-ae  -^1  -ae"^aj2  +  a2e~-^a!l  e-^w2)  (14) 
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Fig.  1.  The  iterative  computation  can  be  modeled 
as  a  first-order  feedback  loop  in  the  iteration 
variable  i. 


Since  H(lo^,cu2)  is  a  separable  filter,  it  is  simple  to  demon¬ 
strate  that  H  (oj^  ,(^2 )  is  stable  if  and  only  if  |a|  <  1. 

However , 


C(to^,w2) 


ae 


—  jui-  .  _  — jaj_  2  -iu>, 

J  1  +  ae  J  2  -a  e  "  1 


e-]w2 


so  that 


(15) 


C(tt,tt)  =  -2a-a^  (16) 

In  particular,  for  a=0.5, 

J  C  ( 7T  ,  TT )  |=  1.25  >  1  (17) 

which  violates  the  convergence  criterion. 

In  general,  if  A(w1,w2)  is  a  trigometric  polynominal 
given  by  equation  (2)  and  C ((0^,10 2)  is  a  trigometric  polynominal 
given  by  equations  (4)  and  (5),  then  H(u)^,u)2)  will  be  continuous 
and  finite  as  long  as  Cfw^,^)  ?  1.  In  this  case,  the  2-D 
inverse  Fourier  transform  can  be  applied  to  H(u>^,u)2)  to  yield 
the  impulse  response  h(n^,n2),  which  will  satisfy 


Z  Z  |h  (n1  ,n2)  |  <  00 


nl  n2 


(18) 


Consequently,  the  condition  C(u1,co2)  /  1  is  a  sufficient  test 
for  BIBO  (bounded  input-bounded  output)  stability  of  H(u^,u>2) 
since  we  are  not  concerned  with  the  region  of  support  of  the 
impulse  response  h(n1,n2).  Clearly,  if  the  convergence 
criterion  |  C  (u^  ,<d2  )  |  <  1  is  satisfied,  then  H(to1,w2)  is  stable 

since  C(w^,oo2)  ?  1.  The  converse,  however,  is  not  true,  as 
the  counter-example  has  shown. 
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D.  TERMINATION  OF  THE  ITERATIVE  COMPUTATION 


In  any  practical  implementation,  the  number  of  iterations 
actually  computed  is  finite.  Using  equation  (9) ,  it  is  straight¬ 
forward  to  show  that  the  output  signal's  spectrum  after  the 
Ith  iteration  is  given  by 


Yi(o)1,o)2)  =  A(a)1,u2)X(a)1,u)2) 


I 

Z  C1  (u>.  ,o)_) 
i=0  i 


=  Yfu^,^)  [l-CI+1(a)1 


f  w2)  ] 


(19) 


We  can  take  the  ratio  Y^  (oj^  ,oj2) /Y  ,  oj2  )  as  one  measure  of 

the  spectral  error  introduced  by  terminating  the  iterative 
computation  after  I  iterations.  The  ratio  is  complex,  in 
general.  Thus,  for  this  ratio  to  lie  close  to  unity,  it  is 
necessary  to  restrict  |  C  ( , u>2 )  1  so  that 


Yj  ( w  r  ^  2  ^ 
Y  (o)1 ,  u>2 ) 


C(w1,oj2) 


1+1 

<  e 


(20) 


for  all  frequencies  (w^,u>2),  where  e  is  a  small,  positive  con¬ 
stant. 

If  we  specify  a  tolerable  degree  of  spectral  error  by 
fixing  e,  we  can  use  this  relation  to  tell  us  how  many  iterations 
will  be  needed  for  a  given  value  of  |c(w1,u2)|.  Conversely,  it 
can  be  used  to  determine  how  jc(u>^,u)2)  |  must  be  restricted,  as 
in  a  filter  design  algorithm,  if  the  number  of  iterations  I  is 
specified  beforehand. 
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In  addition,  £  can  be  used  to  compute  the  space  domain 
difference  between  the  output  signal  at  the  Ith  iteration 
and  the  desired  output  signal.  It  is  straightforward  to  show 
that 

|y(nrn2)  -  yI(n1,n2)  |  <  eY  (2] 


where 


TT 


-TT 


Although  the  bound  (21)  may  not  be  of  practical  use  since  Y 
may  not  be  known  beforehand,  it  nevertheless  demonstrates 
that  the  sequence  of  output  signals  {y^n^n^}  converges 
to  the  desired  output  signal  y(n^,n2)  when  the  convergence 
criterion  (11)  is  satisfied. 

The  effective  frequency  response  when  the  iterative 
computation  is  terminated  after  I  iterations  is  given  by 

I 

Hj  (u>i  ,oj2)  =  A(u)^f^2^  ^  C  ((*)  ^  ,  u>2 )  (23) 

If  we  neglect  the  factor  A(u^,ui2)  for  the  moment,  Hj(w^,u>2) 
has  the  form  of  a  McClellan  transformation  filter  [4,5]  whose 
tap  weights  are  all  equal  to  one.  HI(o)^,(o2)  can  be  realized 
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by  the  structure  shown  in  Figure  2,  which  embodies  the  iterative 
computation  as  the  signal  passes  from  one  stage  of  the  structure 
to  the  next.  Figure  2  serves  to  illustrate  the  fact  that 
HI(a)^,o)2)  represents  the  frequency  response  of  a  2-D  FIR  filter. 
Thus,  the  filter  response  actually  realized  when  the  iterative 
computation  is  terminated  is  an  FIR  approximation  to  the  desired 
rational  frequency  response  H{io^,u)2). 

In  the  McClellan  transformation,  the  kernel  Cfco-i^a^)  is 
usually  associated  with  a  particular  class  of  2-D  FIR  filters 
such  as  those  exhibiting  approximate  circular  symmetry.  The  tap 
weights  are  varied  to  realize  the  mapping  of  a  particular  1-D 
FIR  filter  into  its  2-D  counterpart.  However,  since  the  itera¬ 
tive  computation  has  the  same  form  at  each  iteration,  the 
structure  shown  in  Figure  2  must  be  identical  from  one  stage  to 
the  next,  forcing  all  the  tap  weights  to  be  equal.  The  kernel 
C  (u)^,w2)  itself  is  varied  to  approximately  realize  the  desired 
frequency  response  H  (o^  ,(u2)  . 

Let  us  define  h^n^f^)  to  be  the  effective  2-D  finite- 
extent  impulse  response  implemented  by  I  iterations  of  the 
iterative  computation.  Because  of  the  nature  of  the  iterative 
computation,  the  effective  impulse  response  grows  in  extent  at 
each  iteration.  For  example,  if  we  assume  A(iii^,u)2)=l  and 
c(nlfn2)  has  an  NxN  region  of  support,  then  when  the  iterative 
computation  is  terminated  after  I  iterations,  hj(n^,n2)  will 
have  an  INxIN  region  of  support.  The  number  of  multiplications 
required  to  implement  hj(n^fn2)  by  the  iterative  computation 

however,  is  only  IN^,  in  general.  Although  the  area  of  the  non- 

,  2 
zero  region  of  hI(n1,n2)  goes  up  proportionally  as  I  ,  the  in¬ 
crease  in  the  number  of  multiplications  is  proportional  to  I. 
Since  the  terminated  iterative  implementation  can  be  viewed  as 
a  special  case  of  the  McClellan  transformation  implementation, 
the  more  detailed  results  of  [5]  can  be  applied  to  questions  of 
computation  and  storage. 
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EFFECTS  OF  FINITE  PRECISION  ARITHMETIC 


To  examine  the  effects  of  finite  precision  arithmetic 
on  the  accuracy  of  the  output  signal,  let  us  return  to  the 
consideration  of  the  iterative  computation  as  the  first-order 
feedback  loop  shown  in  Figure  1.  The  true  output  spectrum 
Y(w^,oj2)  will  satisfy  the  implicit  equation 

Y  (o)^ ,  u>2 )  =  A(u)1,(D2)X{a»1,u2)  +  C  (w. ,  w  )  Y  ,  w2 )  (24) 

At  each  iteration,  noise  will  be  injected  into  the  computation. 
We  shall  assume  that  this  noise  occurs  due  to  fixed-point 
round-off  at  the  summation  box  in  Figure  1.  For  a  particular 
ordered  pair  (n1#n2),  the  convolution  of  c(n1,n2)  with 
yi-l^nl'n2^'  rePresented  by  the  operator  C(u>^,a)2),  and  the 
subsequent  addition  to  the  corresponding  point  of  the  signal 
a (nlfn2) *x (n^ ,n2)  can  be  performed  with  double-precision 
arithmetic  and  the  result  rounded  back  to  single  precision. 
Thus,  at  each  iteration,  the  computation  yields 

%  'X, 

Yi((D1,w2)  =  A  (u)^  ,u>2 )  X  (lo^  ,a)2 )  +  C  (a)^,0)2)  Yi_1  (w1,a)2) 

+  ei(ci)1,o)2)  (25) 

'Vj 

where  Y^(o)^,u)2)  is  the  spectrum  of  the  signal  actually  computed 
on  the  ith  iteration  and  e^(a>^,w2)  is  the  spectrum  of  the  com¬ 
putation  noise  introduced  at  the  ith  iteration.  If  we  define 

A  % 

E^{ultu2)  =  Y  (uj^,u>2)-Y^  (co^,u>2) 
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and  subtract  equation  (25)  from  equation  (24)  ,  we  get 

E  i(u1,0)2)  =  C(o)1,(D2)E^_1(aJ1,aj2)  +  e^tu^,^)  (26) 

which  can  be  interpreted  as  a  first-order  difference  equation 
in  the  iteration  variable  i  with  input  (u)^,oj2)  and  output 
Ei(a)^,a)2).  Under  the  assumption  that  the  noise  spectrum 
ei(u)^,w2)  is  uncorrelated  from  iteration  to  iteration,  the 
techniques  used  in  chapter  9  of  [6]  can  be  applied  to 
equation  (26)  to  give  the  expression  for  the  output  noise 
power  under  the  assumption  that  C(u>^,u>2)  satisfies  the  con¬ 
vergence  criterion  (11).  Thus,  if  a^( uj^,a>2)  represents  the 
noise  power  at  (ui^,u>2)  due  to  arithmetic  round-off  and 
aE2  ((i>i  ,oj2)  represents  the  noise  power  in  the  output  spectrum, 
then 


aE  (ti)1,u2)  * 


a^~  (u^  ,u>2) 


1  -  I C  <01^ , u>2 )  I 


(27) 


If  we  further  assume  that  a2  (u^,^)  is  a  constant  and  equal  to 
P  then  the  total  noise  power  (TNP)  in  the  output  is  given  by 


TNP 


4ir‘ 


It 

f  f  - - - *—  dco.dco^ 

JJ  1- 1  C(u)1,uj2>  | 2  1  2 


(28) 


—  IT 
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Intuitively,  we  would  expect  that  if  C ( , aj2 )  satisfies 
the  convergence  criterion  (11)  ,  then  the  effects  of  the 
round-off  noise  e^tio^,^)  introduced  at  the  ith  iteration 
would  gradually  die  out  as  the  iterative  computation  continues. 
Of  course,  additional  round-off  noise  is  introduced  at  each 
iteration,  and  equation  (28)  gives  the  total  noise  power  in 
the  output  signal  at  the  steady  state. 

If  we  define  c. (n, ,n«)  as  the  inverse  Fourier  transform 

l  1  2 

of  C1  (  '  t^ien  equation  (28)  can  be  re-written  using 

Parseval's  theorem. 

*  o 

TNP  =  P  l  E  l  ci(n1,n2)  (29) 

i=0  nx  n2 


When  the  iterative  computation  is  terminated  after  I  iterations, 
the  summation  over  i  in  equation  (29)  extends  only  to  i=l.  This 
result  is  similar  to  that  obtained  in  [5]  for  the  total  noise 
power  of  a  McClellan  transformation  implementation  like  Figure  2. 
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F.  A  SIMPLE  EXAMPLE 


Suppose  we  want  to  implement  a  2-D  digital  filter  whose 
frequency  response  is  approximately  circularly  symmetric 
and  is  approximately  equal  to  a  given  real,  symmetric  1-D 
rational  function  along  each  frequency  axis.  The  McClellan 
transformation  can  be  applied  to  the  1-D  rational  function  to 
generate  a  real,  symmetric  2-D  rational  frequency  response 
which  can  be  implemented  approximately  by  the  iterative 
computation.  For  example,  consider  the  1-D  rational  function 


H  (oj ) 


A(ai) 

I-C(uj) 


where 


4 

A (u)  *  a(0)  +  I  2a  (m)  cos  mu> 
m=l 


2 

C(w)  =  E  2c  (n)  cos  nw 
n=l 


with 


a (0)  *  0.50033 
2a(l)  =  0.81561 
2a (2)  =  0.41543 
2a (3)  =  0.11311 
2a (4 )  =  -0.000037952 


and 


2c(l)  *  -0.0010153 
2c (2)  -  -0.83047 
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This  function  is  real  and  symmetric,  and  it  has  the  shape  shown 
in  Figure  3.  When  the  McClellan  transformation  for  circular 
symmetry  is  applied  to  this  function,  we  get  the  real,  sym¬ 
metric  2-D  rational  function  shown  in  Figure  4.  This  function 
has  the  familiar  form 

H(co1,oj2)  =  A(oj1,o)2)/[1-C  (w1,a)2)  ] 

and  furthermore,  C(cj^,co2)  satisfies  the  convergence  criterion 
(11).  Figure  5  gives  an  indication  of  the  form  of  C(u>^,a>2)  by 
showing  its  value  along  the  oi^  frequency  axis. 

After  20  iterations,  the  surface  plot  of  the  effective 
frequency  response  of  the  iterative  implementation  is  not  dis¬ 
tinguishable  from  Figure  4  because  of  the  resolution  limitations 
of  the  plotting  device.  To  observe  the  convergence  of  Hj.  (u>^,u)2) 
to  the  desired  H^u^,^),  it  is  more  informative  to  plot  the 
values  of  along  either  of  the  two  frequency  axes. 

Figure  6  shows  evaluated  along  the  frequency  axis 

for  1=15,20  and  25  iterations. 

When  comparing  Figures  6  and  3,  it  is  easiest  to  see  the 
difference  between  H(u1,o)2)  and  Hi(o)1,oj2)  at  the  origin,  where 
H(o)1,o)2)  and  |  C  ( co^ ,  u>2 )  |  are  large.  Here,  |C(0,0)|=  0.83,  so 
that  after  15  iterations,  the  relative  error  e  defined  in 
equation  (20)  and  expressed  as  a  percentage  is  approximately 
6.1%.  After  20  iterations,  e  is  approximately  2.4%,  and  after 
25  iterations,  e  is  less  than  1%. 
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Fig.  3.  The  real,  symmetric  1-D  rational  frequency  response  H(u>) 
is  shown  for  0  £  oj/2tt  £0.5.  It  also  represents  the  values  of 
H  shown  in  Figure  4  along  either  of  the  two  frequency 

axes. 


Fig.  5.  The  function  C(uj_,ui2)»  evaluated  along  either  of 
the  frequency  axes,  is  equal  to  C(u>).  It  is  plotted  for 
0  <  oj/2tt  <  0.5. 


G.  DISCUSSION 


The  iterative  implementation  presented  here  has  a  number  of 
interesting  facets.  It  can  be  used  to  implement  a  2-D  rational 
frequency  response  in  an  approximate  manner  by  decomposing  it 
into  a  sequence  of  FIR  filtering  operations.  Thus,  non-causal 
2-D  HR  filters  can  be  realized  approximately  by  this  implemen¬ 
tation  without  performing  the  computationally-intensive  task  of 
spectral  factorization  (or  partial  fraction  expansion) . 

The  rational  frequency  response  to  be  implemented  does  not 
need  to  be  tested  for  BIBO  stability.  Instead,  it  must  satisfy 
a  convergence  criterion  which  is  somewhat  more  restrictive  and 
which  implies  BIBO  stability. 

The  terminated  iterative  implementation  can  be  regarded  as 
a  special  case  of  a  transposed  direct  form  implementation  of  a 
McClellan  transformation  filter  (Figure  2) ,  which  in  turn  is  a 
generalization  of  a  transposed  direct  form  tapped  delay-line 
filter  [5]  .  in  the  terminated  iterative  implementation,  the 
tap  weights  all  turn  out  to  be  equal  to  one.  The  structure  of 
figure  2  is  constrained  to  be  identical  at  each  stage,  since  the 
iterative  computation  has  the  same  form  at  each  iteration. 

Because  of  this  property,  the  terminated  iterative  implementation 
does  not  have  quite  the  generality  of  the  McClellan  transformation 
implementation.  On  the  other  hand,  the  fixed  tapped  weights  do 
not  enter  in  any  design  algorithm  minimizations,  nor  do  they 
require  an  additional  multiplication  of  the  2-D  signal  by  a  constant 
at  each  iteration. 

By  viewing  the  iterative  computation  as  the  first-order 
feedback  loop  shown  in  Figure  1,  we  are  led  to  two  possible 
generalizations  of  the  iterative  implementation.  The  first  is 


to  consider  more  complicated  signal  flow  diagrams  such  as  a 
second-order  feedback  loop.  The  second  is  to  consider  an  input 
sequence  of  2-D  signals  which  may  vary  in  some  way  from  one 
iteration  to  the  next. 

The  iterative  implementation  provides  one  method  of  approx¬ 
imately  implementing  a  real,  symmetric  2-D  rational  transfer 
function  with  a  device  that  can  realize  small-extent  FIR  filter 
kernels.  It  also  provides  an  interesting  theoretical  bridge 
between  2-D  HR  filters  and  2-D  FIR  filters  which  may  be  of  use 
in  conceptualizing  relationships  between  these  two  predominant 
class  of  filters. 
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